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Large transverse momentum distributions of identified particles observed at RHIC are analyzed 
by a relativistic stochastic model in the three dimensional rapidity space. Temperature for inclusive 
reactions is estimated. 
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I. INTRODUCTION 

At RHIC colliding energy of nuclei becomes up to 200 
GeV, and thousands of particles are produced per event. 
To describe such many particle system, a sort of collec- 
tive approach will be useful. One-particle rapidity or 
pseudo-rapidity distributions observed at RHIC are well 
described by the Ornstein-Uhlenbeck process . 

As for the transverse momentum px distribution ob- 
served at RHIC, it has a long tail in the GeV region 
compared with an exponential distribution in px- 

In reference , an empirical formula for large pt dis- 
tributions at polar angle 9 = tt/2, 
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was proposed from the analogy of Landau's hydrodynam- 
ical model. In Eq.Q, E denotes energy of an observed 
particle, Lt is a parameter, and j/t is called the trans- 
verse rapidity. Equation well describes the large px 
distributions for p + p — > ir° + X and p + p — > + X . 
However, it cannot be derived from the hydrodynamical 
model. 

The transverse rapidity is defined in the geodesic cylin- 
drical coordinate system in the three dimensional rapid- 
ity space . The Lorentz invariant phase volume element 
in it is given as 
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= m 2 sinh £ cosh £ dyd^dtj). 
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In the above equation, y denotes the longitudinal rapid- 
ity, and £ denotes the transverse rapidity; 
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where E,pl, Pt and m denote energy, longitudinal mo- 
mentum, transverse momentum, and mass of the ob- 
served particle, respectively, and m-r = yPr + 771-2 ■ ^ 
should be noted that yr coincides with £, only if 9 = tt/2. 

We have proposed the relativistic diffusion model, and 
analyzed large pt distributions for charged particles in 
Au + Au collisions Q . The distribution function of it is 
gaussian-like in radial rapidity, and resemble with Eq. IpLjl 
at 9 = tt/2. 

In section 2, the relativistic diffusion model is briefly 
explained. In section 3, analyses of large pt distribu- 
tions for identified particles, 7r°, 7t~, K~ and p observed 
at RHIC 0, IE @ are made. Temperature is also esti- 
mated from pt distributions. Final section is devoted to 
summary and discussions. 



II. DIFFUSION EQUATION IN THE THREE 
DIMENSIONAL RAPIDITY SPACE 

For simplicity, we consider the diffusion equation with 
radial symmetry in the geodesic polar coordinate system, 
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with initial condition 

^ ^' ^ Att sinh^ p 

In Eq.®, D is a diffusion constant, and p denotes the ra- 
dial rapidity, which is written with energy E, momentum 
p and mass m of observed particle, 
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Inversely, energy and momentum are written respectively 
as 



p\ + Pt = 777 smn P- (5) 



E = m cosh p, |p 

The solution 0, E| of Eq. J5J with the initial condition 
© is given by 
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A physical picture described by Eqs.10) and 10 is as 
follows; After a collision of nuclei, particles are produced 
at the origin of rapidity space expressed by Eq. © . Then 
those particles propagate according to the diffusion equa- 
tion (J2J). In the course of the space time development, en- 
ergy is supplied from the leading particle system (collided 
nuclei) to the produced particle system. Number density 
of particles becomes lower and at some (critical) density, 
interactions among secondary particles cease, and parti- 
cles become free. 

We can analyze transverse momentum (rapidity) dis- 
tributions at fixed polar angle 9, using the equation, 



f(p,t) = C-i— 
sinn p 

a{tf = 2Dt, 



exp 
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with parameters, C and <r(t) 2 , 0]. where transverse mo- 
mentum is given by |px| = msinhpsinfl 



III. ANALYSIS OF p T DISTRIBUTIONS 
OBSERVED AT RHIC 

Transverse momentum (pr) distributions of identified 
particles observed by the PHENIX collaboration [J 0. E| 
are analyzed. The results on px distribution in p + p — > 
7T° + X is shown in Fig. Solid curve shown in Fig.l 
is drawn by the use of Eq.©, parameters of which are 
estimated with the least mean square method, and are 
shown in Table 1. The results on pr distributions in 
Au + Au — > 7T° + X are shown in Fig. [3] and Table 2. 
Observed px distributions on 7r° both in p+p and Au+Au 
collisions are well described by Eq.©. 
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FIG. 1: p T distribution for p + p -> n° + X at y = H 

The results on Au + Au — > tt~ + X are shown in Fig. [31 
and Table 3. As is seen from Fig. 3 and Table 3, fitting 



TABLE I: Estimated parameters on pr distributions in p + 
p -> 7T° + X at y = at </s = 200 GeV Q 
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a{tf 


Xmin/md.f. 


1388.0 ± 121.5 


0.602L 0.004 
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FIG. 2: p T distribution for Au + Au -> tt° + X at y = 



of the theoretical curve to the observed px distribution 
is gradually improved as the centrality cut increases, in 
other words, the average number of participant nucleons 
decreases. 

The results on Au + Au — » K ~~ + X are shown in Fig.0] 
and Table 4, and those in Au + Au — > p + X are shown 
in Fig. Eland Table 5. As can be seen from Figs. 4 and 
5, and Tables 4 and 5, fitting of the theoretical curve 
to the observed pr distribution are better than those in 
Au + Au — > 7r~ + X, and become much better as the 
centrality cut increases. 



TABLE II: Estimated parameters on pr distributions in Au+ 
Au -> ty° + X at y = at yfs = 200 GeV @ 



centrality 


C 


a(tf 


Xmin /fl.d.f 


00-10% 


16088 ± 2003 


0.574 ± 0.005 


27.3/16 


10-20% 


1161.5 ± 1390.1 


0.580 ± 0.005 


20.0/15 


20-30% 


7345.9 ± 922.7 


0.586 ± 0.005 


14.9/13 


30-40% 


4378.8 ± 569.3 


0.593 ± 0.005 


15.7/13 


40-50% 


2559.8 ± 345.1 


0.600 ± 0.006 


9.4/12 


50-60% 


1425.3 ± 194.3 


0.599 ± 0.007 


12.4/12 


60-70% 


527.1 ± 78.9 


0.617 ± 0.008 


8.2/10 


70-80% 


241.2 ± 39.6 


0.616 ± 0.009 


7.2/9 


80-92% 


123.0 ± 24.7 


0.611 ± 0.011 


5.0/8 
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FIG. 3: p T distribution for Au + Au — > 7r + X at y = @ 



FIG. 4: p T distribution for Au + Au — > if + X at y = 



TABLE III: Estimated parameters on pr distributions in 
Au + Au -> tv~ + X at y = at y/s = 200 GeV @ 



Cl± L 1 CX11 by 


C 


a(tf 


xLn/n.d.f 


00-05% 


2539.6 ± 88.8 


0.701± 0.007 


522.8/26 


05-10% 


2092.7 ± 73.7 


0.710± 0.008 


447.4/26 


10-15% 


1692.0 ± 59.9 


0.718± 0.008 


390.7/26 


15-20% 


1391.1 ± 49.9 


0.724± 0.008 


319.4/26 


20-30% 


1052.8 ± 37.5 


0.727± 0.008 


278.5/26 


30-40% 


667.0 ± 24.4 


0.736± 0.009 


144.8/26 


40-50% 


409.8 ± 15.4 


0.736± 0.009 


144.8/26 


50-60% 


241.2 ± 9.3 


0.731± 0.010 


114.4/26 


60-70% 


125.4 ± 5.1 


0.724± 0.010 


94.3/26 


70-80% 


57.4 ± 2.5 


0.719± 0.012 


60.4/26 


80-92% 


29.5 ± 1.4 


0.708± 0.013 


42.0/26 


TABLE IV: 


Estimated parameters on pr distributions in 


Au + Au — » 


K ~ + X at y = 


at y/s = 200 GeV 


a 


centrality 


C 


a(tf 


xLn/n.d.f 


00-05% 


44.75 ± 2.81 


0.456± 0.019 


9.45/14 


05-10% 


37.05 ± 2.33 


0.457± 0.019 


8.67/14 


10-15% 


30.82 ± 1.93 


0.453± 0.018 


10.77/14 


15-20% 


25.69 ± 1.62 


0.45 1± 0.018 


8.52/14 


20-30% 


19.17 ± 1.21 


0.450± 0.018 


7.07/14 


30-40% 


12.92 ± 0.82 


0.437± 0.017 


8.37/14 


40-50% 


7.86 ± 0.50 


0.431± 0.017 


6.57/14 


50-60% 


4.55 ± 0.31 


0.420± 0.018 


3.85/14 


60-70% 


2.27 ± 0.16 


0.409± 0.019 


4.44/14 


70-80% 


1.01 ± 0.08 


0.396± 0.021 


3.60/14 


80-92% 


0.50 ± 0.04 


0.386± 0.025 


7.44/14 
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FIG. 5: pr distribution for Au + Au -^p + X&ty = 



In order to estimate the temperature for inclusive re- 
actions from our analysis, we consider the approximate 
expression for Eq.@) in the small |px| region. When 
p « 1, sinhp ~ p. Then, Eq.© reduces to 



f(p,t) = Cexp 
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2a(ty 



(7) 



Equation Q should coincide with the Maxwell distribu- 
tion, 
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TABLE V: Estimated parameters on pr distributions in Au+ 
Au -> p + X at y = at y/s = 200 GeV @ 



centrality 




c 








a(tf 




xLn/n.d.f 


00-05% 


6.48 


± 





29 





296± 





008 


100.1/20 


05-10% 


5.49 


± 





25 





296± 





008 


84.31/20 


10-15% 


4.57 


± 





21 





297± 





008 


75.89/20 


15-20% 


3.85 


± 





18 





296± 





008 


65.83/20 


20-30% 


3.04 


± 





14 





287± 





008 


64.19/20 


30-40% 


2.10 


± 





10 





281± 





008 


39.75/20 


40-50% 


1.30 


± 





07 





277± 





009 


21.33/20 


50-60% 


0.79 


± 





05 





266± 





010 


10.40/20 


60-70% 


0.40 


± 





03 





254± 





012 


2.23/20 


70-80% 


0.18 


± 





02 





239± 





016 


6.29/20 


80-92% 


0.10 


± 





01 





225± 





019 


3.18/20 



The estimated temperature kT from the px distribu- 
tion for 7T° is about 80 Mev, and that for ir~ is about 100 
MeV. These temperatures are almost independent from 
the centrality cut. On the other hand, the temperature 
for K~ distributions gradually decreases from 230 MeV 
to 190 MeV, as the centrality cut increases from 5-10% 
to 80-92%. That for p distributions decreases from 280 
MeV to 210 MeV as the centrality cut increases. 

Our analysis suggests that K~ and p produced at the 
lowest centrality cut, 0-5%, would keep somewhat earlier 
memory than 7r°-meson or 7r~-meson. It is very inter- 
esting to know whether identical particle correlations of 
if -mesons, protons and so forth depend on the centrality 
cut or not. 



IV. SUMMARY AND DISCUSSIONS 



400 — , — , — | — , — , — , — | — , — , — , — | — , — , — , — | — , — ,— 

Centrality cut dependence of kT 

350 - 

a kT (pbar) 

300 - ° "Oil 

.... ° «!*) 



□ nun □ 




150 - 

100^2OOOO O O O O O O 



50 - 

t_. — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i— 
20 40 60 80 100 

central ity cut (X) 



FIG. 6: Centrality cut dependence of temperature kT esti- 
mated from pt distributions in Au + Au collisions 



Then we have an identity, 

kT = ma(t) 2 . (8) 

From Eq.||SJ, we can estimate the temperature kT of 
inclusive reactions for observed particle with mass m. 
The results are shown in Fig. EJ We use m^o = 135, 
m^- = 140, m K - = 494, and m p = 938 MeV. 



In order to analyze large pr distributions of charged 
particles observed at RHIC, a stochastic process in the 
three dimensional rapidity space is introduced. The so- 
lution is gaussian-like in radial rapidity. 

Transverse momentum distributions for 7r° and 7r~ at 
y = at yjl = 200 GeV observed by the PHENIX col- 
laboration are analyzed. Observed px distributions for 
7T° in pp and Au + Au collisions are well described by the 
formula JHJ), and the result of fitting in Au + Au collisions 
becomes much better as the centrality cut increases. 

The result on transverse momentum distributions for 
7r~ in Au + Au collisions are not good as those for ir°. 
However, as the centrality cut increases, the result be- 
comes much improved. This tendency would suggest 
that the initial condition J3J) is simpler for the heavy 
ion collision process. At the lower centrality cut, col- 
lisions among projectile and target particles occur much 
more times compared with those at higher centrality cut. 
Therefore, secondary particles produced at each collision 
may have a certain (collective) velocity distribution wider 
than the delta function © even in the initial stage. In 
order to include the effect, we should change the initial 
condition ©, or consider the diffusion equation without 
radial symmetry in the three dimensional rapidity space. 
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